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Abstract. This paper reviews several analytic tools for the field of topological 
insulators, developed with the aid of non-commutative calculus and geometry. 
The set of tools includes bulk topological invariants defined directly in the 
thermodynamic limit and in the presence of disorder, whose robustness is shown 
to have non-trivial physical consequences for the bulk states. The set of tools 
also includes a general relation between the current of an observable and its edge 
index, relation that can be used to investigate the robustness of the edge states 
against disorder. The paper focuses on the motivations behind creating such tools 
and on how to use them. 
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1. Introduction 

The topological insulating state represents an entirely novel manifestation of the 
ordinary matter, characterized by the existence of metallic edge states despite the 
presence of strong disorder or strong deformations of the materials. We will argue 
here that the topological insulating state is also characterized by the existence of 
bulk states that resist localization. These are intrinsic properties of the topological 
insulators and no external magnetic held is need to trigger this peculiar behavior [T]. 

Fervent theoretical activity around the topological insulators resulted in a variety 
of characterizations, some based on very original ideas [21 El HI El [HI El El [SI EH [HI HH 
H31 H31 H3 [111 EZl HB]. For example, the Z 2 invariant introduced by Kane and Mele 
[2] was new to the physics and mathematics communities. Same can be said about 
the Spin-Chern number introduced by Sheng et al [10] or about the recent results on 
3D topological insulators. We also want to mention that several topological invariants 
appeared long before the held of topological insulators was defined, in the context of 
3 He-A like topologically ordered systems [T5] . 

Most of the analytic works were completed within the framework of translational 
invariant bulk systems and homogeneous edges, where one can build topological 
invariants on the Brillouin torus and the edge states can be computed explicitely. 
The topological properties in the presence of disorder were mostly investigated by 
computational means [101 Ell E21 E%1 El HH1 HH EU- The present contribution to this 
special issue on "Topological Insulators" gives a short overview of the author's recent 
works [25l [26l [271 El] , which provide a set of analytic tools to diagnose the effect of 
disorder on the bulk and on the edge states of the topological insulators. 

The topic of the present paper is how to define topological invariants directly in 
the thermodynamic limit and in the presence of disorder and what kind of information 
can be extracted from their invariant properties. The discussion excludes technical 
arguments, which can be found in the author's published work, but includes a careful 
presentation of the assumptions and the limits in which the results are valid. The 
connection to other works, especially to the existing numerical studies, is emphasized. 
The paper also pays special attention to the motivation behind these works and how 
to use the emerging analytic tools. The author also shares some of his current personal 
beliefs that crystalized during this work. 

2. Tools for the bulk states 

We start the discussion from the classical Chern invariant. For any smooth family of 
projectors {Pe}oeM indexed by 9, who lives on a two dimensional smooth and closed 
manifold M, one can compute the Berry curvature 

dF = (2TTi)- 1 Ti{Pg[de 1 Pe,de 2 Pe}}d9 1 A d6 2 (1) 

and then integrate the curvature over M to obtain the Chern number, an integer that 
is invariant to smooth deformations of the family of projectors. We want to emphasize 
from the beginning that the Chern number can be defined only if Pg is at least first 
order differentiable in 6, for all 8 of the manifold M. The Chern number can change 
its value if Pg seize to be a global smooth function of 0. 

The last two observations are quite relevant because many invariants for 
topological insulators were linked to the Chern numbers of families of projectors 
built through various procedures, without any check on the global smoothness of 
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with respect to k (the Bloch wave- vector) . Without such check, the constructions 
remain formal and incomplete. The issue should not be treated lightly since proving 
global smoothness of various objects on the Brillouin torus represents an outstanding 
problem in the band theory. Physicists as well as mathematicians bumped into it 
when trying to extend Kohn's results on the analytic structure of Bloch functions 
from one to higher dimensions (29] [30l EU 1221 EI]. What usually happens is the 
following: the objects are easily seen to be locally smooth on k and, starting from a 
small region of the Brillouin torus, the objects can be smoothly extrapolated to larger 
and larger regions. But when it is to cover the entire torus, point singularities in the 
first derivatives can develop that cannot be removed by any gauge transformation. 
Proving the absence of such singularities is almost always a formidable task [30] . 

The construction of the spin-Chern number given in Ref. [28J includes a check of 
the global smoothness of the projectors involved in the construction. Let us repeat the 
construction of Ref. [5H] , this time emphasizing this check. We will use the concrete 
model of electrons in graphene [Ml [5] : 

H = -t \ m , cr )( n , '\ 

(mn) ,a 

+i^so [s • (dfe m x d n h)]<T,a'\m,a)(n,a'\ , 2 ) 

{(mn)) .era' 

+i\r [z • (s X d m n)]a,<T>\n,o-)(n,o-'\. 

{mn) .aa' 

Here, s is the spin operator, m and n denote the sites of the honeycomb lattice and a 
and a' the electron spin degrees of freedom, taking the values ±1. The Hamiltonian 
acts on the Hilbert space TL spanned by the orthonormal basis \n, a). The simple 
angular brackets in Eq.[3]denote the nearest neighbors and the double angular brackets 
denote the second nearest neighbours. Inside the second sum, k represents the unique 
common nearest-neighbor of m and n. Also, d mn is the displacement from the site m 
to the site n, with the units chosen so that the distance between nearest neighbors is 
unity. The electrons are considered non-interacting. The three terms in Eq. [2] are the 
usual nearest neighbor hopping term, the intrinsic spin-orbit coupling preserving the 
lattice symmetries and the Rashba potential induced by the substrate supporting the 
graphene sheet. We assume that the parameters in the Hamiltonian arc chosen so that 
we are in the Spin-Hall zone of the phase diagram. At some point in the paper, we will 
view the honeycomb lattice as a triangular lattice of unit cells containing two C atoms. 
When we do so, we will use the notation \i,a.,o~) for the elementary quantum states, 
where i represents the position of a unit cell in the triangular lattice and a = 1, 2 is 
a label for the two atoms inside each unit cell. 

There are two sites per unit cell and two spin states per site, therefore the model 
has 4 states per unit cell. For the beginning, we will stay with the translational 
invariant case, where we can use the Bloch fibration, i.e., the unitary transformation 
from the Hilbert space of the infinite sample into a continuum direct sum of 4 
dimensional complex spaces: 

U:H^® keT C\ (C/*)(fc) = ^E J -e-* j *(j), (3) 

where k lives on the Brillouin torus T = [0, 27r] x [0, 2w]. U transforms the original 
Hamiltonian in a direct sum of Bloch Hamiltonians: UHqU^ 1 = (B^eT Ho(k), which 
display two upper and two lower bands separated by an insulating gap. Let P denote 
the spectral projector onto the states below the insulating gap of Hq. Under the 
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Bloch fibration, this projector becomes: UPU~ 1 =Q) keT P(k) and the Chern number 
associated with the family of projectors P(k) is zero, as it will generically be for any 
time reversal invariant band model. According to Ref. [35], the fiber bundle of the 
occupied states is topologically equivalent to the trivial fiber bundle TxC 2 . 

It is clear that, in order to obtain non-trivial Chern numbers, we need to split the 
projector P(fe)=P_ (fe)©P + (fe) and hope that P±(k) have non-trivial Chern numbers. 
Due to the time-reversal invariance and the half-integer value of the spin, one can 
always split P(fe) like this and there are in fact an infinite way of doing the splitting 
at each k. This observation was used, for example, in Refs. [T7] and [T5] to give 
alternative methods of calculus for the Z2 invariant. But the big question is: Among 
this infinite ways of splitting, is there at least one that leads to P±(k) that are globally 
smooth over the entire Brillouin torus? 

It appears to us that, for the model of Eq. [21 the only way to construct such 
splitting is to use the s z operator (s z \i,a, a) — \o~\i, a, <r)). The key idea is to 
use the spectral properties of the operator Ps z P. Indeed, after the Bloch fibration, 
UPs z PU~ x becomes fce7 - P(k)s z P(k) and we can diagonalize each of the operators 
P(k)s z P(k). If the Rashba term doesn't exceed a threshold value, the spectrum of 
P(k)s z P{k) consists of two isolated eigenvalues, positioned symmetrically and away 
from the origin, for all the fc's of the Brillouin torus. The spectral projectors P±(k) 
onto the positive/negative eigenvalues are smooth of k and can be used to achieve the 
decomposition P(k) = P-(fc) © P + (k). At this point we can can define the Chern 
numbers C± for the families of projectors P±(k) and define the spin-Chern number 
as C S =\(C + — CL). For the model of Eq. [21 C s =l (assuming A50 positive). 

The smoothness of P±{k) follows from a direct computation for the particular 
model of Eq.[21 but a more general, model independent, proof was achieved in Ref. [2B] 
based on the fact that global smoothness of P±(k) implies exponential decay of the 
real space kernel P±(i, j) with the separation \i— j\, and viceversa, exponential decay 
implies global smoothness. The exponential decay of P±(i,j) can be established using 
traditional methods [55] , 

Besides allowing for an explicit check of the global smoothness, the above 
construction also allows one to write the spin-Chern number without using the 
Brillouin torus, something that is absolutely necessary when dealing with disorder. 
Indeed, the Chern numbers C± can be computed in real space using the following 
formula, which is completely equivalent to the k space formula: 

C± = 2vrz tr{P± [[fi,P±], [x 2 ,P±]]}, (4) 

where "tr" refers to the trace over the four quantum states of the first unit cell, [, ] 
denotes the usual commutator and x is the position operator: x\i,a, a) = i\i,a,o~). 
P± are the spectral projectors onto the positive and negative parts of the spectrum 
of Pa z P, respectively. This formula can be used to define the Chern numbers for the 
disordered case. Indeed, if we add a random potential V u to Hq, one can compute 
the projector P w onto the occupied states of H UJ =Ho+V L j and then investigate the 
eigenvalue spectrum of P w <t z Pj. For \r not exceeding a threshold value, one will 
find that the spectrum segregates into two isolated islands, one above and one below 
zero. Therefore, one can define the spectral projectors P^ onto the upper/lower 
islands of spectrum of Pu& z P u and compute C± using Eq. [4j A fundamental result in 
noncommutativc geometry |36] states the following: 

Proposition 1. Let H u be a family of projectors depending on the random variable 
oj. If the matrix element decays sufficiently fast with the separation \i — j\ 
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(see the following discussion) and if IT/s satisfy u^II^m" 1 = n triW , with u n being the 
lattice translation by an arbitrary n and t n a flow in the space of the random variable, 
then: 

C = 2m (tr{n w [[ii,n w ], [£ 2 ,IL]]}) W , (( ) u = disorder average) (5) 

is an integer that is invariant to smooth deformations of H w 's as long as they remain 
localized. 

Ref. [28 showed that, if P^ is exponentially localized and the gap of P w s z P w 
remains open, then P^- are exponentially localized. The gap of PuSzPu was shown to 
remain open as long as is smaller than a threshold value and the disorder is not very 
strong. Moreover, since the spin operator commutes with translations, the relation 
UnP^u^ 1 — Pt u is automatically satisfied, therefore P^ satisfy the conditions of 
Proposition 1 and consequently C± and C s are topologically invariant integers as long 
as P u is exponentially localized and the gap of PuiS z Pu, remains open. 



2.1. Physical consequences 

The spin-Chern number can be used for topological classification but it will be equally 
important to extract non-trivial physical consequences from its remarkable properties, 
in particular, to say something about the localization of the bulk states in the presence 
of disorder. For the model defined by Eq. [2l we can establish the following fact: a 
non-zero spin-Chern number implies existence of energy regions where the localization 
length of the bulk states diverges, at least for small Xr. In other words, there are bulk 
states that resist localization in the presence of disorder. 

Indeed, let us imagine a numerical experiment in which we lower the Fermi level 
continuously, from the mid-gap all the way to -co, and try to repeat the construction 
of the spin-Chern number, this time starting from P w (E F )a z P w (E F ) \P u (Ep) = the 
projector onto the energy spectrum below Ep\. When the Fermi level enters the bulk 
energy spectrum, the spectrum of P^(Ef)o' z P u1 (Ef) can rapidly spread within the 
entire [— h, |] interval, engulfing the origin. Even though, according to the general 
results of Ref. [36], the Chern numbers for projectors P^(Ep) remain well defined 
and invariant of Ef as long as the spectrum of P u {Ef)o~ z Puj{Ef) remains localized 
near the origin, more precisely as long as: 

(Ti{k P±{E f )x 2 P±{E f )M) u < oo, (6) 
where 7To is the projector onto the four states of the first unit cell. This condition says 
that the kernel 

P^E F (0,a,a;i,a',a') = (0,a,o-\P±(E F )\i,a',a') (7) 
decays fast enough so that: 

1/2 



^(Ef 



Ew 2 EE i^ BF (o,w,«>')i 5 



(8) 



is finite. The quantity inside the square root is precisely the left side of Eq.[Hl A ± (Ep) 
can be viewed as (polarized-) localization lengths. 

Our real interest is in the true localization length: 

A(E F ) - [{Tt^P^Ef^P^Ef)^}) J 1/2 . (9) 
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The connection between the true and the polarized localization lengths is given by the 
following relation: 

K{E F f = K+{E F f +A-{E F ) 2 

+ (Tr{w P+(Ep)x?p-(Ep)wo}) u (10) 

+ (Tr{^ Pj(£ F )x 2 P+(i; F h }) w • 

But the last two lines are null when A/j=0, and they remain small relative to the terms 
appearing in the first line, for small A^. 

Now as we move the Fermi level to — oc, the Cern numbers C± become zero, 
therefore they must jump at some point from ±1 to zero. This can happen only 
if A ± (Ep) diverge as Ep crosses certain energy regions. According to Eq. [TO] this 
implies that the true localization length A(Ep) also diverges when Ep crosses these 
energy regions, at least for small \r. 

2.2. Discussion 

We would like first to relate our findings to the numerical studies of Refs. [TOl [24] . In 
the following, if not otherwise specified, the Fermi level is considered in the middle 
of the bulk gap. The first observation is that the original construction in Ref. [28] 
considered same type of disorder as in these references, a white noise with amplitude 
between [-W/2.W/2] (to keep the notation of Refs. PH El]]). For the model of 
Eq. 1, by direct computation in the clean limit, we have verified that the gap of 
PuS z P u remains open when \r is increased, all the way until the gap of Hq closes. 
This tells us that C s is a robust topological invariant, for the entire domain of QSH 
phase, in total agreement with the numerical observations of [TU]. In Ref. [55], the 
amplitude W was assumed small so that a gap remains open in the spectrum of H^, 
so this work established the robustness of C s only in the presence of relatively small 
disorder. The study by Sheng and her collaborators went beyond the small disorder 
case and established that C s remains well defined even if the insulating gap is filled 
with localized spectrum. 

Switching to the second numerical study [24] . we would like to specify that the 
discussion of Ref. [25J was limited to the range W < 3 of Ref. [24} . In this case, the 
exponential localization of P w and of P^ can be explicitly proven; therefore, according 
to our arguments, one should observe bulk states that resist localization, in line with 
the numerical findings of Ref. [24] , Beyond W = 3, the spin-Chern number remains 
a topological invariant only if the Fermi level can be placed in a region of localized 
states, as it was already discussed above. Referring to Fig. 1 of Ref. [23], this is the 
case for panels (a-1), (b-1) and (c-1), which correspond to W = 5 and different values 
of A_r, but not for panels (a-2), (b-2) and (c-2) which correspond to W=7. C s can be 
defined again for panels (a-3), (b-3) and (c-3), correspoinding to W=8, but this time 
it takes a trivial value. It is interesting to notice that the evolution of the de-localized 
spectral regions with the increase of W seen in these numerical experiments can be 
actually predicted using the following arguments. As long as the Fermi level can be 
placed in a region of localized spectrum, the spin-Chern number remain non-trivial, 
therefore we should observe de-localized states below and above Ep. It is known [37] ■ 
however, that for very large values of W, all the states localize, therefore at some 
point there should be no alternative, but to place Ep in a de-localized region of the 
spectrum. But Ep can be always placed in a region of localized spectrum unless the 
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de-localized spectral regions below and above Ep move towards each other and at 
some point "scissor" the Fermi level. This is precisely what was observed in Ref. |24J. 

Secondly, we would like to discuss some predictions. The Z 2 topological invariant 
can be computed as C s (mod2). According to our arguments, there should be bulk 
states that resist localization not only for odd but also for even spin-Chern numbers, 
that is, for trivial Z2 invariants. This, for example, must be the case for the model of 
Eq. [2] with spin |, which has |C S |=2. Therefore, it will be very interesting to see the 
calculations of Ref. [21] repeated for this system. 

Another observation is that the spin-Chern number can be defined even in the 
absence of time-reversal symmetry. Ref. [24] presented a calculation that included 
a Zeeman term in the x direction, of strength h x . The results for relatively large 
magnetic field indicate that the bulk states are entirely localized (although the 
graph corresponding to the weaker disorder is not very conclusive in our opinion). 
Unfortunately, Ref. [51] does not discuss what happens at weaker magnetic fields, 
i.e. when the time reversal symmetry is only weakly broken. In this regime, the last 
two lines of Eq. [10] remain small and a divergence of A ± (Ep) implies a divergence of 
A(Ep). In other words, we predict that, for small values of h x , there are bulk states 
that resist localization in the presence of disorder. 

Our final remark for the section is that the present construction of the spin-Chern 
number can be repeated for other operators, let us call them w, instead of s z . If the 
eigenvalue spectrum of PwP consists of two or more isolated islands, one can compute 
the Chern numbers of the spectral projectors for each island and generate w-Chern 
numbers. A non-zero ui-Chern number is a signal of topological properties, leading to 
a possible zi-Hall insulator. We hope that this observation will lead to more efficient 
search and discovery paths for new classes of topological insulators. 

3. Tools for the edge states 

In 2004, the following result appeared [38 . Consider the Landau Hamiltonian 
describing free electrons in the semiplane z=0 and x>0 (therefore with an edge along 
y axis), moving in a periodic potential, a uniform magnetic field along the z axis (no 
rational flux required) and a random potential V^. Then the following is true: 



where p(e) is any arbitrary distribution with support in a bulk energy gap of H u and 
such that J p(e)de=l, the function F is defined by F(e)=J^° p(s)ds and x+( x iV) 1S 
the characteristic function for the semi-plane y>0. The trace tro is taken over the 
states with support in the plane y—0. 

This is a remarkable result since it contains the first explicit proof that the 
edge states of the Landau Hamiltonian resist localization in the presence of disorder. 
Indeed, let us read again Eq. [11] The left hand side is the total charge current of the 
edge states inside a bulk energy gap of H^, weighted with the distribution p. This 
current will be zero if the edge states localize. Therefore, if one can show that the right 
hand side is non-zero, he has an explicit proof that the edge states resist localization. 
But the right hand side is the index of an operator (to be discussed shortly) , which can 
be computed explicitly and for the Landau Hamiltonian it turns out to be a non-zero 
integer number. One can see then how Eq. [Tl] becomes a very effective tool for the 
study of the edge physics. 




(11) 
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The above result is relevant for Landau Hamiltonians and for Chern insulators 
[26] , but unfortunately is not relevant for time-reversal invariant topological insulators 
because their edge charge current is zero. Therefore, we need to look at the currents 
of other observables (spin, etc.). This was the main motivation behind the work of 
Ref. [25]. Using non-commutative calculus, the author established a general result 
stating when and why is the current of a general observable X quantized. The result 
can be summarized as follow. Assume: 

a) Existence of a family of self-adjoint Hamiltonians H u on 7i, with u> S fl. 

b) Existence of a 1-parameter unitary group u a such that UaXu^ 1 = X + a. 

c) Existence of an ergodic flow t a on Vt such that UaH^u^ 1 = H tabJ for all ui 6 CI. 

Proposition 2. If ir + denotes the projector onto the positive eigenvalue spectrum of 
X, p{e) is a statistical distribution of the quantum states such that J p(e)ale=l and 
F(€)=$™ p{e)de, then: [below, X u (t) = e itH - Xe~ itH -] 

(tr |p(JT„)^^ }} = ^Indexj^e-^^V} , (12) 

provided the kernel of e~ 2mF ( H "> has certain localization properties. 

The index of an operator is equal to the difference between the number of zero modes 
of that operator and the number of zero modes of its conjugate. It is an integer 
number, invariant to continuous deformations of the operator. The index on the right 
hand side of Eq.rT2]is now called the edge index of observable X, and can be explicitly 
computed by deforming the Hamiltonian into simple, exactly solvable models. During 
such deformations, one has to make sure that the technical conditions are not violated, 
that is, the localization of e - 27 ™ F (^). 

Eq [12] can be a very useful tool for the analysis of the edge spectrum. Note 
that the left side of Eq. rT2] is zero if the spectrum of is completely localized. 
Therefore, a non-zero value of the right side of Eq . [T2l assures that has absolutely 
continuum spectrum, even in the presence of disorder. It becomes clear then that 
Eq. [12] establishes a very general method to demonstrate the robustness of the edge 
spectrum against disorder, which amounts to a search of an observable X for which 
the index of the right hand side of Eq. rT2]is non-zero. 

For Quantum Spin-Hall insulators, this observable was defined in Ref. [27J. 
The construction of the observable X was based on the fact that for spin i 
particles and time- reversal invariant Hamiltonians, the Hilbert space can be split 
into two sectors H=TC^®H+ such that each of the sectors are left invariant by the 
Hamiltonian and r7i.± = H.^ (t = time- reversal operation) . The observable X is simply 
X=P + yP + —P_yP-, where P± are the projectors onto the sectors H± and y is the 
coordinate along the edge. For the model of Eq. [2] with an edge and random potential, 
the edge index of this observable is equal to spin-Chern number. In other words, 
using the above tools, we now have an explicit proof that the edge states of Quantum 
Spin-Hall insulators resist localization when disorder is present. 

The original and now most popular argument [34] for the robustness of the edge 
states is based on the observation that the scattering between the right and left moving 
edge bands cannot occur in time-reversal invariant insulators. In author's opinion, this 
argument is not enough. The reason is because, although we are concerned with the 
quasi-one dimensional edge states, we are still dealing with a 2D systems. In such 
systems, a band can scatter into an infinite number of analytically continued bands 
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as detailed in Refs. [3U] and [35] ■ In other words, if one takes an energy e in the 
bulk insulating gap and e n ^ denotes the band energies of the model with an edge, 
the equation e=e n ^ has an infinite number of solutions, thought many of them will 
occur at complex k. Nevertheless, the incoming edge state ip ln will scatter into these 
complex k states, which are the analytic continuation of the bulk states, and time- 
reversal invariance cannot prevent this. For these reasons, it appears to the author 
that the edge bands can backscatter through the bulk. 

4. Outlook and concluding remarks 

We hope that the present paper makes more clear the motivation behind the somewhat 
abstract tools developed with the aid of non-commutative calculus and geometry. 
They allow one to define robust topological invariants directly in the thermodynamic 
limit and in the presence of disorder. Their robustness can have direct physical 
implications and the non-commutative methods can lead to precise statements about 
the localization of the bulk and the edge states. The list of known topological insulators 
is expected to increase in the following years and the generalization of the spin-Chern 
number to the w-Chern number could aid the search and discovery of new classes 
of topological insulators. The edges of these new topological insulators could cary 
currents of new observables, therefore the generality of the tools developed for the 
edge physics may become more useful in the future. 

One of the interesting and important things left unsolved is how to define the Z2 
bulk topological number without involving the Brillouin torus and how to construct 
its non-commutative theory? What are the physical implications of its robustness 
against disorder for the bulk states? In the author's view, these are the most pressing 
questions in the field of topological insulators. 

As one can see, the discussion of the paper did not touch the three dimensional 
case at all. Unfortunately, the applications of non-commutative geometry to 3D 
systems are very scarce. It is author's hope that the discovery of the 3D topological 
insulators and the recent and exciting theoretical work on these systems will renew 
the interest of non-commutative geometers and new applications will come soon. 
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